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SUMMARY 


An  experimental  and  analytical  study  has  been  made  of  some  features 
of  the  turbulent  heat  diffusion  Behind  a  line  heated  -wire  stretched 
perpendicular  to  a  flowing  isotropic  turbulence.  The  mean  temperature 
distributions  have  been  measured  with  systematic  variations  in  wind 
speedy  size  of  turbulence -producing  grid^  and  downstream  location  of 
heat  source.  The  nature  of  the  temperature  fluctuation  field  has  been 
studied. 

A  comparison  of  Lagranglan  and  Eulerian  analyses  for  diffusion  in 
a  nondecaying  turbulence  yields  an  expression  for  turbulent -heat -transfer 
coefficient  in  terms  of  turbulence  velocity  and  a  Lagranglan  "scale." 

The  ratio  of  Eulerian  to  Lagranglan  microscale  has  been  determined 
theoretically  by  generalization  of  a  result  of  Heisenberg  and^  with 
arbitrary  constants  taken  from  independent  sources^  shows  rough  agree¬ 
ment  with  experimental  results. 

A  convenient  form  has  been  deduced  for  the  c:);*iterion  of  inter¬ 
changeability  of  instantaneous  space  and  time  derivatives  in  a  flowing 
turbulence. 


INTRODUCTION 


One  of  the  most  striking  aspects  of  turbulent  motion  in  fluids  is 
its  dispersive  property.  This  "convective  diffusion/Villustrated  by 
the  general  statistical  tendency  of  (noncontiguous)  fluid  elements  to 
get  farther  apart  with  increasing  time^  was  probably  first  observed  long 
before  the  era  of  analytical  fluid  mechanics.  An  analytical  start  on 
this  problem  was  not  raade^  however^  until  the  now-classic  work  by 
Taylor  in  1921  on  diffusion  by  continuous  movements  (reference  l) .  Not 
only  did  this  paper  lay  a  groundwork  for  the  study  of  turbulent  diffusion 
but  it  also  represented  a  forward  step  in  the  ideas  essential  to  develop¬ 
ment  of  a  general  statistical  theory  of  turbulence^  a  field  which  had' 
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scarcely  progressed  since  Reynolds’  original  formulation  of  the  equations 
of  motion  for  a  flow  in  which  mean  and  fluctuating  parts  could  he 
distinguished. 

The  diffusive  action  of  a  turbulent  flow  may  manifest  itself  in 
various  ways^  depending  upon  the  initial  and/or  boundary  conditions  and 
upon  the  interests  of  the  observer.  The  following  possible  measures  of 
the  diffusive  powers  are  neither  exhaustive  nor  mutually  independent: 

(1)  The  average  rate  of  dispersion  of  particles  from  a  fixed  source 

(2)  The  average  rate  of  increase  of  spacing  between  different 

particles 

(3)  The  average  rate  of  transport  of  particle  concentration  under 

a  given  mean  concentration  gradient 

(4)  The  average  rate  of  increase  of  the  length  of  a  fluid  line 

(5)  The  average  rate  of  increase  of  the  area  of  a  fluid  surface 

The  word  "particles”  means  simply  indelibly  tagged  fluid  elements^  much 
smaller  than  the  smallest  length  associated  with  the  turbulence. 

The  present  report  is  concerned  primarily  with  measure  (l).  The 
measurements  have  all  been  made  in  the  thermal  wake  of  a  long  thin  heated 
wire  mounted  perpendicular  to  an  isotropic  turbulent  air  flow  and 
producing  no  turbulent  wake.  Here  the  tagging  is  thermal^  and  the  degree 
of  indelibility  (negligibility  of  molecular  diffusion)  is  one  of  the 
matters  to  be  investigated. 

The  diffusive  property  (for  scalars)  of  a  turbulent  flow  is  appar¬ 
ently  a  secondary  characteristic  at  least  in  the  sense  that  it  need  not 
explicitly  enter  the  dynamical  problem.  The  diffusion  may  be  regarded 
as  a  kinematic  phenomenon^  to  be  deduced  from  the  dynamical  solution  to 
the  problem  if  and  when  the  latter  is  obtained.  Thus  the  objective  of 
research  on  turbulent  diffusion  may  be  to  seek  a  connection  between  the 
diffusive  and  the  dynamical  statistical  variables,  even  before  the 
complete  dynamical  theory  is  available. 

Measure  (3)  is  usually  termed  the  "turbulent  transport"  or  "transfer" 
problem.  Although  of  extensive  practical  importance,  it  has  not  yet  been 
subjected  to  genuine  theoretical  study. 

Most  of  the  semiempirical  "theories"  of  turbulent  transport,  for 
both  scalar  and  vector  properties,  employ  an  Eulerian  formulation  of  the 
basic  equations,  and  up  to  now  they  have  been  \mable  to  relate  the 
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turbulent  transport  correlation  to  other  statistical  functions  describing  ^ 
the  flow.  Taylor  (reference  l)  showed  that  in  the  simple  case  of  a 
homogeneous  field  of  isotropic  turbulence^  and  even  in  a  decaying 
isotropic  turbulence  (reference  2)^  a  Lagrangian  formulation  of  the 
transport  (i.e._,  diffusion;  the  terms  will  be  used  interchangeably) 
problem  leads  to  some  important  results.^ 

Up  to  the  present  time  little  theoretical  or  experimental  work  has 
'been  done  to  find  relations^  if  any^  between  the  Lagrangian  statistical 
measures  of  a  turbulent  field  with  its  Euler ian  statistical  measures. 

Since  turbulence  dynamics  seems  best  handled  in  the  latter  terms  and 
turbulent  diffusion  in  the  former^  it  is  evident  that  such  a  connection 
is  important.  Hence^  one  of  the  purposes  of  the  present  experiments  has 
been  to  compare  the  magnitudes  of  some  of  these  quantities  under  varia¬ 
tions  in  the  turbulent  field.  For  example^  the  postulates  of  Taylor 
and  Heisenberg  on  a  relation  between  Lagrangian  and  Eulerian  microscales 
can  be  examined  and^  in  corrected  form_,  compared  with  experiment. 

The  turbulent  diffusion  from  a  fixed  line  source  can  be  set  up 
analytically  as  an  ordinary  (Eulerian)  ^’heat-transfer"  problem,  permitting 
a  start  to  be  made  in  relating  measures  (l)  and  (3)  of  the  diffusive 
power  of  a  turbulent  flow,  under  certain  simplifying  assumptions. 

Measures  (4)  and  (especially)  (5)  may  well  be  classed  as  character¬ 
istic  of  the  "turbulent-mixing"  problem  rather  than  diffusion  in  the 
common  connotation. 

Experimental  work  on  diffusion  from  a  fixed  local  source  in  a 
turbulent  flow  has  been  meager.  In  isotropic  turbulence,  there  have 
been  the  measurements  of  Schubauer  (reference  4),  Simmons  (reported  by 
Taylor  in  reference  2),  Dupuis  (reported  by  Kampe  de  Feriet  in  refer¬ 
ence  5)^  Frenkiel  (reference  6),  and  Collis  (reference  7)*  Of  these, 
only  the  data  of  Simmons  and  Collis  are  extensive  enough  to  permit 
confident  computation  of  the  Lagrangian  correlation  function.  In 
turbulent  shear  flow,  Skramstad  and  Schubauer  (reference  8),  Dryden 
(reference  9)^  aiid  the  present  authors  (reference  lO)  have  measured 
distributions  close  to  a  source;  Kalinske  and  Pien  (reference  11)  and 
Van  Driest  (reference  12)  have  made  measurements  somewhat  farther 
downstream. 

None  of  these  studies  was  repeated  with  a  systematic  variation  of 
the  properties  of  the  turbulence.  In  spite  of  the  poor  precision 


lln  his  tensorlal  generalization  of  Taylor's  work  on  case  (l), 
Batchelor  (reference  3)  has  chosen  to  call  this  an  "Eulerian"  analysis, 
describing  only  case  (2)  as  "Lagrangian."  In  keeping  with  previously 
accepted  nomenclature,  both  cases  are  Lagrangian,  (l)  involving  a  single 
particle,  and  (2)  involving  a  pair.  In  fact,  case  (2)  might  be  termed 
a  mixed  (Eulerian  and  Lagrangian)  problem. 
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Inherent  in  this  type  of  measurement^  it  was  hoped  that  such  an  approach 

would  at  least  show  up  some  general  trends  in  the  relations  between 

Eulerian  and  Lagrangian  variables.  w 

This  investigation  has  been  conducted  at  the  Aeronautics  Department 
of  the  John  Hopkins  University  under  the  sponsorship  and  with  the 
financial  assistance  of  the  National  Advisory  Committee  for  Aeronautics. 

The  authors  would  like  to  acknowledge  the  assistance  of  Messrs.  Alan 
Kistler^  George  Stierhoff,  and  Allen  Gates  and  Miss  Patricia  O’Brien^ 
as  well  as  the  helpful  criticism  of  Dr.  Francis  H.  Clauser  and  Dr.  C.  C. 

Lin. 


SYMBOIS 


X 

Xo 

Ax  =  X  -  Xq 

y 

z 

r 

U 

V 
u 

V 

w 


c 


0 


grid  mesh  size 

distance  downstream  from  grid 
location  of  heating  wire 

distance  in  direction  of  measured  diffusion 
distance  in  direction  of  heating  wire 
scalar  distance  between  two  points 
mean  velocity  in  x-direction 
mean  velocity  in  y-direction 

instantaneous  velocity  fluctuation  in  x-direction 
instantaneous  velocity  fluctuation  in  y-direction 
instantaneous  velocity  fluctuation  in  z -direction 

root -mean- square  molecular  velocity 

instantaneous  temperature  (measured  above  ambient  room 
temperature) 
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0 

00 

=/P 


,5 

mean  temperature 

maximum  mean  temperatiire  at  a  cross  section^  a  function 
of  /SK 

instantaneous  temperature  difference 
temperature 'difference  of  rectangular  heat  pulse 


j  width  of  rectangular  heat  pulse 

s  average  on-center  spacing  of  pulses 

rt 

Tl(to,t)  =  V(|)  dl 

lo 

t  time 


T  time  difference 

^*(^)  Eulerian  velocity  correlation  coefficient  (notation  of 

Von  Karman  and  Howarth)  * 

Lagrangian  correlation  coefficient  for  nondecaying 
turbulence 


Lagrangian  correlation  as  a  function  of  •  t)  for 
nondecaying  and  decaying  t'urbulence 


distance  traveled  in  y-direction  by  a  fluid  particle 
distance  traveled  in  x-direction  by  a  fluid  particle 

root -mean- square  displacement  of  a  fluid  particle  in 
y-direction 
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root -mean- square  displacement  of  a  molecule 

/  / 

Eulerian  microscale  of  turbulence  X  = 


1 


Lagrangian  microscale  of  turbulence  for  nondecaying 


turbulence  (Xy  = 


1/2^ 


-Rv"(0) 


Lagrangian  microscale  of  turbulence  for  nondeca;;^ -^ag  and 


and  decaying  turbulence  =  Lr  q ) 


1/2^ 


Poo 


Lagrangian  scale 


Rv(t )  dr 


Eulerian  scale  IL  =  f(r)  dr 

V  .  0 


Lagrangian  scale  for  nondecaying  and  decaying  turbulence 

dTij 

mean  free  path  of  a  molecule 
mixing  length 

specific  heat  per  unit  volume  at  constant  pressure 
density 

kinematic  viscosity 
thermal  conductivity 
static  pressure 

dimensionless  empirical  constant 
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poo 


H  =  2pCpU 


©  dy 


H*  = 


kip 


H 


2pcp  f(|)  d| 

Jo 


turbulent-heat-transfer  coefficient 


Mt  = 


■n  BTi(n) 


T  = 

p(^) 

Rl 


probability  density  of  fluctuating  temperature 
turbulence  Reynolds  number  based  on  Eulerian  scale 


R 


X 


n 

max 

min 


turbulence  Reynolds  number  based  on  Eulerian  microscale 


mean  value  or  ensemble  average 
maximum^  used  as  subscript 
minimum^  used  as  subscript 
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EQUIPMEOT  AKD  PROCEDURE 
Aerodynamic  Equipment 


The  wind  tunnel  (fig.  1)  is  an  open-return  NPL  type  tunnel  with 
a  2-  by  2-foot  working  section  and  a  free-stream  turbulence  level  of 
v’/u  =  0.06  percent  and  u*/u  =  0.05  percent  at  a  mean  velocity  of 
26  feet  per  second.  The  turbulence -producing  grids  were  as  follows: 


Designation 

Type 

Mesh  size 
(in.) 

Rod 

diameter 

(in.) 

Solidity 

l“in.grld 
1/2-in.  grid 
l/4-ln.  grid 

Biplane,  wood 
Woven,  metal 
Woven,  metal 

1.00 

^.50  by  0.53 

.25 

0.25 

.120 

.060 

0.1^37 

.41 

.42 

^The  0.50“ln.  mesh  was  set  in  direction  of  measured  diffusion. 


They  were  mounted  in  txirn  at  the  upstream  end  of  the  working  section. 

The  heat  source  was  an  0.008-inch-diameter  platinum  wire  stretched 
vertically  across  the  tunnel  at  various  distances  from  the  grid.  It  was 
heated  by  direct  current  to  temperatures  between  500°  and  700°  C^  with 
the  latter  figure  only  at  the  highest  operating  velocity  of  38.0  feet 
per  second.  The  wire  Reynolds  numbers  at  this  condition  and  at  the  two 
other  velocities  were  as  follows: 


Mean  speed 
(fps) 

Wire  temperatizre 
(°C) 

Reynolds  number  based  on  - 

Air  temperature 

Wire  temperature 

8.5 

500 

35 

6.4 

25.6 

,  500 

105 

19 

38.0 

700 

157 

19 

A  preliminary  investigation  was  made  without  grids  to  Insure  that 
these  operating  conditions  did  not  generate  a  vortex  street  downstream 
of  the  heated  wire. 

With  the  grids  in  place,  the  mean  momentum  wake  became  practically 
undetectable  with  total -head  tube  and  manometer  at  distances  greater 
than  1  or  2  inches  downstream  of  the  heating  wire. 
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Measuring  Equipment  and  Procedures 

The  mean-temperature  distributions  vere  measured  with  a  Chromel- 
Alinael  thermocouple  and  a  Leeds  and  Northrup  type  K-2  potentiometer. 

The  cold  Junction  was  kept  outside  of  the  wind  tunnel. 

The  shadowgraph  technique  was  used  to  photograph  the  laminar  thermal 
wake  close  to  the  soxirce  with  no  grid  in  the  wind  tunnel.  This  informa¬ 
tion  was  applied  to  the  problem  of  "correcting"  the  thermal  wake  in 
turbulent  flow  for  the  effects  of  molecular  diffusion  and  of  finite 
source  size.  Figure  2(a)  is  a  shadowgraph  of  the  wire  wake  with  no  grid 
in  the  tunnel}  figxire  2(b)  is  a  typical  time  exposure  with  grid-produced 
turbulence.  A  resistance -thermometer  traverse  of  the  laminar  wake  in  a 
flow  of  very  small  tiirbulence  showed  that  the  temperature  profile  had 
already  become  very  nearly  Gaussian  at  a  distance  of  1  inch  (125  wire 
diameters)  downstream.  The  white  lines  on  the  sides  of  the  dark  wake 
shadow  in  figure  2(a)  correspond  to  the  minimums  in  the  second  derivative 
of  the  density  profile.  For  small  temperature  differences  these  coincide 
with  the  maximums  in  the  second  derivative  of  the  temperature  profile. 
Although  the  temperature  differences  are  not  small  in  the  Immediate 
vicinity  of  the  wire^  this  condition  is  reasonably  well  satisfied  at 
relatively  small  values  of  Ax  as  evidenced  by  the  parabolic  spread  of 
this  laminar  wake. 

The  standard  deviation  of  the  wake  could  then  be  computed  from  the 
spacing  of  these  two  bright  lines,  with  the  assumption  of  a  Gaussian 
distribution.  Since  the  closest  points  of  traverse  in  the  turbulent 
cases  were  l/2  inch  (63  wire  diameters)  from  the  heat  source,  this  was 
probably  a  reasonable  assumption. 

As  pointed  out  by  Taylor  (reference  2),  the  molecular  and  the 
turbulent  diffusive  phenomena  are  statistically  independent,  so  that  the 
squares  of  the  standard  deviations  due  to  these  two  effects  are  additive. 
Hence  the  wake  spread  due  to  turbulence  alone,  from  a  true  line  source, 
was  obtained  by  subtracting  the  square  of  the  standard  deviation  of  the 
laminar  wake  (computed  from  the  shadowgraph)  from  the  square  of  the 
standard  deviation  of  the  total  wake  (from  thermocouple  traverses  with 
grid  turbulence  present)  at  all  stations.  This  difference  was  the 
square  of  the  standard  deviation  of  the  desired  phenomena.  All  wake- 
spread  data  presented  in  the  next  section  have  been  corrected  in  this 
fashion. 

Parenthetically,  it  should  be  remarked  that  the  laminar  wake  in 
figure  2(a)  spreads  parabolically  within  the  limits  of  precision,  from 
at. least  1  inch  on,  so  that  the  effects  of  density  differences  on  the 
flow  phenomena  must  have  been  negligible  for  this  investigation. 
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The  transverse  tiorhulence  levels  /u  behind  the  grids  were 
obtained  from  the  initial  rate  of  spread  of  the  mean  thermal  wake 
(method  due  to  Schubauer  (reference  ^4-))  after  the  effects  of  molecular 
spread  and  finite  source  had  been  removed.  The  resulting  levels  were 
somewhat  higher  than  those  obtained  with  a  hot-wire  anemometer  but  were 
used  because  of  their  consistency  with  the  rest  of  the  measured  dif¬ 
fusion  curve. 

Free-stream  velocity  fluctuations  (without  grids)  and  the  wake 
temperature  fluctuations  (with  grids)  were  measured  with  the  hot-wire 
anemometry  equipment  described  in  reference  10.  The  wires  were 
0.00025-inch  platinum  etched  from  Wollaston;  the  compensated  response 
of  the  system  was  flat  within  l2  percent  over  a  frequency  range  from 
3  to  12^000  cycles  per  second. 

Oscillograms  of  the  temperature  fluctuations  were  recorded  by 
photographing  a  blue  oscilloscope  tube  with  fast  35-railli^eter  film 
in  a  General  Radio  type  65I-AE  camera. 

Probability  densities  of  the  temperature  fluctuations  at  fixed 
points  in  the  mean  thermal  wake  were  determined  from  photodensitometer 
traverses  of  time-exposure  photographs  of  a  short-persistence  ( 0.001  sec 
blue  oscilloscope  tube  with  the  temperature  fluctuations  on  one  pair  of 
plates  and  a  30^000-cycle-per-second  sweep  on  the  opposite  plates.  The 
technique  is  essentially  that  used  by  Simmons  and  Salter  (reference  13). 


EXPERIMENTAL  RESULTS 
Mean  Thermal  Wake 


Complete  mean- temper atirre  wakes  behind  a  line  source  of  heat  were 
measured  for  10  different  conditions.  Arranged  to  indicate  the 
systematic  variation  of  one  parameter  at  a  time^  these  conditions  were 
as  follows: 


l-ln.  grid;  —  =  43.4j  U,  fps,  . 

8.5 

25.6 

38.0 

U  =  25.6  fpsi  ^  =  43.4;  grid,  In.,  . 

1 

1/2 

1/4 

U  =  25.6  ^  =  86.1;  grid.  In., . 

1 

1/2 

1/4 

U  =  25.6  fps;  1/2-in.  grid;  xq/M  . 

43.4 

86.1 

172.3 

U  =  25.6  fps;  l/4-ln.  grid;  Xq/M  . 

43.4 

86.1 

172.3 

Here  U  is  the  mean  velocity^  M  is  the  grid  mesh  size^  and  Xq  is 
the  heat-source  location  measured  from  the  grid.  Since  some  individual 
cases  enter  as  elements  in  two  sequences^  the  total  number  of  elements 
is  more  than  10. 

Two  of  the  many  mean-temperatirre  traverses  in  the  y-direction 
(perpendicular  to  mean  flow  and  to  source  line)  are  shown  in  figure  3 
to  give  an  idea  of  the  amount  of  experimental  scatter.  The  upper 
traverse  was  the  worst  of  the  lot^  even  showing  an  apparent  skewness 
which  was  not  borne  out  by  the  investigation  as  a  whole.  The  lower 
traverse  is  more  nearly  typical  of  the  measured  temperature  distribu¬ 
tions  from  which  the  standard  deviations  of  the  mean  thermal  wake  were 
computed.  By  comparison  with  the  reference  curve,  it  is  seen  to  be 
essentially  Gaussian.  This  was  the  case  for  temperature  profiles  at 
all  stations.  Since  the  virtually  Gaussian  character  of  such  a  wake 
has  already  been  established  by  several  of  the  earlier  publications, 
there  seemed  to  be  no  point  in  reproducing  here  all  of  the  large  number 
of  traverses  measured. 

The  mean- thermal -wake  spread  for  the  10  different  configurations 
studied  is  given  in  figures  4  to  7  as  plots  of  corrected  standard 
deviation  Y*  against  distance  from  the  heat  source  Z^:^.  Each  point 
in  these  figures  corresponds  to  a  complete  transverse  temperature 
traverse. 

In  order  to  have  values  of  transverse  turbulence  level  v*/u 
consistent  with  the  thermal-wake  behavior,  these  values  were  determined 
from  the  initial  angle  of  spread  of  the  corrected  wake  standard  devia¬ 
tion  (references  2  and  4)  instead  of  from  direct  hot-wire  anemometer 
measurements.  The  results  are  plotted  in  figure  8.  Since  these 
represent  an  insufficient  number  of  points  per  grid  to  permit  the 
drawing  of  reliable  curves,  some  simplifying  assumptions  were  made  based 
upon  the  results  of  several  more-detailed  turbulence-decay  investigations. 
(See  references  l4  to  20.  The  pertinent  results  of  references  l4  to  l6 
are  summarized  in  reference  17*)  In  the  light  of  these  papers,  it  was 
assumed  that  the  decay  curves  had  a  common  apparent  origin,  and  that 
this  was  obtainable  by  drawing  the  best  straight  line  for  all  the 
available  points,  independent  of  the  differing  wind  speeds  and  mesh 
sizes.  In  computing  any  individual  Lagrangian  correlation  function 
from  the  corresponding  wake  history,  the  turbulence  decay  rate  was 
assumed  to  be  given  by  the  line  drawn  through  this  common  origin  and 
the  specific  turbulence  value  giving  the  measured  Initial  spread  angle 
for  this  wake.  This  is,  of  course,  a  very  rough  procedure,  but  the 
experimental  scatter  in  this  whole  method  of  determining  Lagrangian 
correlation  functions  is  so  great  that  a  more  extensive  study  of  decay 
(including  the  resolution  of  inconsistencies  between  wake  method  and 
hot-wire  method)  seemed  unwarranted  at  this  time. 
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Temperature  Fluctuations 

Distributions  of  temperature-fluctuation  level  'S' /s  in  the  thermal 
wake  have  been  measured  by  using  the  hot-wire  anemometer  as  a  resistance 
thermometer^  that  is^  at  a  current  low  enough  to  render  the  sensitivity 
to  velocity  fluctuations  negligible  compared  with  the  sensitivity  to 
temperature  fluctuations  (reference  21).  A  representative  distribution 
of  -6  V0  x-direction  is  given  in  figure  9*  Typical  transverse 

distributions  are  given  in  figures  10(a)  and  10(b).  It  is  clear  that 
the  temperature-fluctuation  intensity  changes  very  little  with  increasing 
values  of  A  rough  explanation  of  the  very  high  values  of  -iS*/® 

(compared  with  the  concomitant  turbulence  levels  for  example)  in  terms 
of  the  highly  intermittent  structure  of  the  thermal  wake  has  been  given 
in  reference  10  and  will  be  discussed  in  more  detail  later  in  this 
report.  This  intermitt ency  is  shown  very  clearly  in  figure  11^  a 
series  of  temperature  oscillograms  recorded  at  two  different  positions 
across  the  thermal  wake  for  a  fixed  value  of  Ax  and  at  two  different 
values  of  Ax  with  y  =  0.  The  one-sided  and  pulse  character  of  the 
instantaneous  temperature  at  a  fixed  point  in  space  is  also  demonstrated 
by  its  probability  density. 


THEORETICAL  CONSIDERATIONS 
Homogeneous  Steady  Turbulence  at  Rest 


For  a  nondecaying  incompressible  turbulence  with  no  mean  motion^ 
Taylor  (reference  l)  was  followed  in  getting  an  expression  for  the  mean 
time  rate  of  diffusion  in  the  y-directlon  (say)  from  a  fixed  soirrce  as 
measured  by  the  second  moment  of  the  probability  density  of  the  diffu¬ 
sion^  that  is^  the  mean-square  particle  displacement  Y^(t): 


dt  dt 


=  2v(t) 


The  bar  denotes  ensemble  average. 

Taking  v(t)  inside  the  integral^  interchanging  the  processes  of 
integration  and  averaging,  and  introducing  the  Lagrangian  (auto) 
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correlation  coefficient  (following  the  fluid  particle),  where 
T  =  t  -  tj^, 


v(t)v(t  -  t) 

there  results 


y2  =  2v2 


t 

0  ^ 


PT 


R„(t)  dr  dT 


(1) 


This  is  Taylor’s  form.  Integration  by  parts  yields  a  form  like  that  in 
the  work  of  Kampe^  de  Feriet  (reference  : 


(t  -  t)Rv(t)  dr  (2) 


In  this  Lagrangian  analysis^  v(t)  is  the  velocity  of  a  fluid 
particle  in  the  y-direction  at  time  tj  v(t  -  x)  is  the  velocity  of 
the  same  particle  at  time  t  -  x.  Corresponding  expressions  can  he 
written  for  the  rate  of  diffusion  in  any  direction. 

Diffusion  from  an  infinite  line  source^  the  case  to  he  discussed 
here,  is  a  two-dimensional  problem  in  the  mean,  and,  in  addition  to 
equation  (2), 


(t  -  x)R^(x)  dx  (3) 


A  tensorial  generalization  of  these  concepts  has  been  given  by 
Batchelor  (reference  3)  dealing  with  the  behavior  of  XiXj(t)  where 
and  Xj  are  any  two  of  the  orthogonal  displacements  of  the  particle  at 
time  t. 


It  should  be  noted  that  this  analysis  gives  no  information  on  the 
shape  of  the  probability  density  of  Y(t)  or  of  X(t) .  In  fact  there 
still  exists  no  theory  for  these.  However,  experiments  in  flowing 
turbulence  (the  case  to  be  considered  next)  show  Gaussian  density, 
within  the  experimental  precision,  for  Y(t)  at  all  values  of  t.  . 
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For  this  stationary  random  process,  Taylor  (reference  l)  introduced 
the  concept  of  the  Lagrangian  "scale," 

^00 

Ll  =  R^(t)  dT  (4) 

0 

These  have  the  dimensions  of  time  and  are  characteristic  constants  of 
the  system. 

In  his  later  work  on  the  (Eulerian)  dynamics  (reference  2),  Taylor 
had  occasion  to  introduce  another  measure  of  the  correlation  function, 
which  he  called  the  "microscale."  Applying  the  same  geometrical  concept 
to  the  present  function,  the  Lagrangian  microscale 


is  simply  the  r-intercept  of  the  vertex-osculating  parabola  of  the  even 
function  Rv('T")-  The  kinematic  significance  is  clearly  shown  by  a 
series  expansion  of  v(t  +  t)  in  Ry(T): 
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From  equations  (6)  and  (5)^  the  Lagrangian  microscale  for  v(t)  is 


or 


X 


2  _ 
V  ” 


(7) 


(7a) 


Turbulence  in  a  Flowing  Medium 

The  dictates  of  both  practical  interest  and  experimental  feasibility 
require  analysis  of  the  diffusion  when  there  is  a  mean  velocity  U 
relative  to  the  source.  Since  the  diffusion  phenomenon  is  linear^  the 
probability  density  (mean-concentration  distribution  of  tagged  particles 
in  the  wake)  is  simply  proportional  to  the  superimposed  probability 
densities  of  a  continuous  line  of  sources  moving  with  the  mean  velocity 
TI  with  their  time  (and  space)  origin  at  the  actual  fixed  source.  This 

is  illustrated  in  figure  12(a)  for  ^  ■  <  1.  The  circles  (corresponding 

to  isotropy)  are  the  standard  deviations  of  the  dispersions  that  would 
occur  from  moving  sources.  The ^envelope  of  these  circles  gives  a 
measure  of  the  mean  wake.  It  is  obvious  that  in  general  the  functional 
form  of  the  mean- concentration  distribution  along  a  line  ^  =  Constant 
will  not  be  the  same  as  the  functional  form  of  the  same  quantity  for 
the  individual  source  at  time  t^  that  is^  at  position  x  =  Ut. 

However^  it  is  extremely  likely  that  if  «  1  (fig.  12(b)) 

the  mean  concentration  along  a  line  Ax  =  Constant  becomes  very  nearly 
of  the  same  functional  form  as  that  for  the  source  at  x  =  Ut.  For  a 
zero-correlation  Gaussian  density^  the  equivalence  is  easily 
demonstrable.^ 


^The  most  general  mathematical  restrictions  under  which  the  super¬ 
position  of  a  line  of  identical  densities  will  yield  a  "cross-section^’ 
density  of  the  same  form  has  not  been  studied  here.  It  is  obvious  that 
the  condition  of  statistical  independence  is  sufficient.  The  fact  that 
the  density  of  each  of  the  velocity  components  in  isotropic  turbulence 
has  been  found  to  be  Gaussian  within  the  experimental  precision  seems 
to  show  that  the  equivalence  under  discussion  is  at  least  a  good 
approximation. 
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The  condition  «  1  will  always  occur  at  large  enoixgh  values 

of  Ax  when  ^  <  1.  This  follows  from  the  asymptotic  parabolic 
behavior,  Y*(t)  ~  \ft  (reference  1),  of  the  diffusive  process.  It  will 

V  ^ 

occur  for  all* values  of  Z\x  when  —  «  1*  Therefore,  in  this  case  a 

U 

simple  approximate  space-time  transformation  in  the  mean  is  permissible^ 
and  the  t-variation  in  Taylor's  theory  of  diffusion  by  continuous  move¬ 
ments  becomes  a  variation  of  x/U.  (it  must  be  emphasized  that  the 
foregoing  discussion  does  not  apply  directly  to  the  possibility  of 
applying  a  space-time  transformation  to  the  instantaneous  turbulence 
variables.  This  latter  question  will  be  discussed  later.) 

V ' 

In  the  present  measurements  ^  therefore,  the 

Ax-variation  of  diffusion  gives  an  approximate  measure  of  the  Lagrangian 
correlation  coefficient  (in  time).  Equation  (2)  can  be  rewritten  as; 


"'(i)  ‘  ^  ?  r  •  HI) 


Decaying  Isotropic  Turbulence 

When  the  turbulence  is  decaying  in  time  (similar  to  space  in  the 

flowing  turbulence)  v^  is  no  longer  constant,  and  the  analysis  cannot 
be  carried  out  as  far  as  equation  (l) .  The  same  approach  stops  with 


v'(t)  v*(t  -  t)  dT 


where  the  prime  denotes  root-mean-square  value,  and 


(9) 


R  =  V(t)v(t  -  t) 
t  t-T  v'(t)v'(t  -  t) 


There  is  no  a  priori  reason  to  believe  that  t^t-T  a  function  of  r 
alone,  as  in  the  nondecaying  turbulent  flow. 


NACA  TN  2710 


17 


At  this  point  Taylor  (reference  2)  invokes  the  empirical  fact  that 
over  a  wide  range  of  mean_velocities  (all  of  which  give  essentially  the 
same  distribution  of  v'/u  in  x  behind  a  grid)  the  thermal  wake 
behind  a  line  heat  source  at  fixed  x  appears  to  be  unchanged  in  form, 
within  the  experimental  error.  This  is  consistent  with  dependence  of 
the  diffusive  process  upon  a  variable  of  type 


Therefore,  Taylor  has  postulated  the  unique  dependence  of 

the  variable  t).  With  this  postulate  and  the  space-time  transformation 
valid  for  small  turbulence  level,  he  arrives  at 


or 


2 


(13) 
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A  physical  significance  of  the  length  r|  is  underscored  by  the 
limiting  form  of  equation  (I3)  as  t] — >0  and  - >1.  Then^ 

Y'  =  T] 

so  that  r\  is  a  measure  of  the  lateral  diffusion  that  would  occur  if 
the  lateral  velocity  fluctuation  following  a  particle  v(t)  remained 
perfectly  correlated  but  decreased  in  magnitude  according  to  the  decay 
rate  of  the  turbulence  level* 

In  postulating  Taylor  was  apparently  comparing  only 

diffusive  processes  in  turbulence  fields  with  identical  — (x) .  Of  more 

general  interest  is  the  comparison  of  diffusion  in  fields  with  differing 
turbulence-level  distributions.  Although  such  a  generalized  application 
of  his  postulate  is  doubtless  not  too  well  applicable^  it  is  conceivable 
it  might  have  approximate  success  in  the  more  general  comparison.  For 
both  convenience  and  lack  of  any  obviously  superior  alternative^  his 
suggestion  is  therefore  applied  in  computing  the  results  of  the  measure¬ 
ments  reported  here. 3 

With  the  T)-postulate_,  the  R^^  Lagrangian  correlation  function  can 
be  obtained  from  measurements  of  as  a  function  of  Ax: 


=  1 

2  ^  dx2 


(1*^) 


^After  this  work  was  completed.  Dr.  Batchelor  suggested  an  alternative 
approximate  approach:  In  order  to  construct  a  stationaiy  random  function 
out  of  the  nonstationary  v(t)  one  first  normalizes  the  dependent 
variable  with  its  root-mean- square  value.  (This  has  been  automatically 
accomplished  by  use  of  the  correlation  coefficient.)  If  the  decaying 
quantity  is  assumed  to  maintain  complete  similarity  during  decay,  all 
characteristic  times  (e.g.,  Lagrangian  time  scale  and  microscale)  vary 
in  the  same  way  with  t  and  the  new  independent  variable  is  constructed 
by  dividing  t  by  this  t-variation.  Unfortunately,  this  variation  is 
unknown  a  priori,  so  it  would  be  necessary  to  assume  further  that  the 
Lagrangian  scales  are  directly  proportional  to  the  Eulerian  scales,  about 
which  there  is  previous  experimental  information.  His  work  has  now  been 
published;  see  reference  22. 
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A  scale  and  a  microscale  can  also  be  defined  for 


(15) 


(16) 


With  nondecaying  turbulence,  no  one-parameter  true  (time)  Lagrangian 
correlation  function  exists,  and  the  r| -formulation  is  much  more  con¬ 
venient.  A  further  significance  of  this  variable  will  appear  in  the 
comparison  of  Eulerian  and  Lagrangian  treatments  of  diffusion  from  a 
line  source  in  flowing  turbulence. 


Accelerations  in  Decaying  Itirbulence 

A  sdries  expansion  of  v(t  +  t)  for  decaying  turbulence  will  show 
something  about  the  initial  behavior  of  the  true  (time)  Lagrangian 
correlation  function  and  will  indicate  an  experimental  method  for 
examining  a  hypothesis  of  Taylor  on  the  interchangeability  of  instanta¬ 
neous  time  and  space  derivatives  when  the  turbulence  level  is  low 
(reference  23) . 

Write  the  Lagrangian  correlation  coefficient 


% 


t) 


v(t)v(t  -  t) 
v'(t)v*(t  -  t) 


(17) 


Substitute 


v(t  -  t)  =  v(t) 


into  numerator  and  denominator,  and  restrict  the  analysis  to  small  values 
of  t: 


(18) 


20 


NACA  TN  2710 


Divide  numerator  and  denominator  by  y  expand  the  square  root  in 
the  numerator,  and  keep  terms  in 


Rjt,T)  ~  1  - 


/dv'  \  ^  /dv\ 
ydt  \dt/ 


^  _1  2v^(t) 


For  negligible  decay  rate  this  reduces  to  equation  (6). 

Equation  (19)  shows  that  a  Lagrangian  microscale  defined  by 


(19) 


lim 


x2(t)  T— 


1  -  Rv(t,T) 


L  J 


T=0 


(20) 


is  expressilDle  as 


Xv  (t)  2v^ 


(if  -  (ff)'_ 


(21) 


Introduction  of  Taylor's  ri-postulate  transforms  equation  (l9)  to 


Rti(^)  1  - 


(22) 


after  the  additional  approximation  that  v'(t])  ~  v'(0)  when  q  is  very 
small. 

This  gives  a  new  expression  for  Lagrangian  microscale 


1  ^  1 
X,l2  =  2^ 


< 
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Equation  (23)  is  ^in  contradiction  to  equation  (17)  of  part  IV  of 
reference  2.  In  that  work  Taylor  has  apparently  assumed  that  v*  is  a 
nondecaying  function  of  t\  I  However  it  certainly  is  decaying_,  even  in 
terms,  of  this  distorted  coordinate^  and  the  dy*/dT]~term  must  he 
included. 


Since  and  v’(x)  can  he  determined  experimentally  from  the 

mean  thermal  wake  behind  a  line  heat  source^  equation  (23)  permits 


determination  of 


"Stokes"  acceleration 


which  is  simply  related  to  the  mean- square 

This  quantity  is  of  particular 


(i)  - 


interest  for  the  possibility  of  an  instantaneous  space-time  transforma¬ 
tion  at  low  turbulence  levels.  This  was  first  proposed  by  Taylor  (refer¬ 
ence  23)  and  has  since  been  used  veiy  widely^  especially  to  get  approxi¬ 
mate  values  of  partial  derivatives  with  respect  to  x  (the  mean-flow 
direction)  by  measurements  of  time  partial  derivatives. 


The  total  (or  Stoke s_)_  derivative  of  v(x^yjZ,t)  in  a  turbulent 
flow  with  mean  velocity  U  along  the  x-direction  is 


dv  _  ^v 
dt  ~  St 


(U  +  u) 


Sx 


+  w 


Sz 


(24) 


Taylor's  hypothesis  amotmts  to  the  statement  that 


§t 


or 


with 


Sv 

St 


ox 


(24a) 


U2 


1/2 


«  1 


(25) 
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In  detail^  equation  (25)  is 


u  -u . 

1  J 


5v  bv 
bx±  Sxj 


2u- 


dv  bv 
dt  bx^ 


1/2 


«  1 


¥ 


(25a) 


In  the  absence  of  information  on  the  algebraic  sign  of  the  triple 
correlation  term^  it  is  sufficient  to  require  the  two  conditions^ 


IS 


m 


1/2 


«  1 


bv  bv 
^  J  ^xi  ^ X  ^ 


u-u 


1/2 


«  1 


(26) 


But 


(iJ 


can  be  determined  from  measurements  of  Y’(x)} 


/N  \  2  2 

\S^/  ”  ^  where  X  is  the  Eulerian  microscale  (reference  2)j  and 

upper  bounds^  in  terms  of  measurable  functions^  can  be  set  on 


Uj^U4  with  the  use  of  Schwarz  inequalities.  Thus,  an  experimental 

check  of  the  requirements  in  equation  (26)  is  to  be  made  in  the  section 
entitled  "Computation  of  Results." 


^Lin  has  discussed  the  validity  of  Taylor ^s  hypothesis  using  a 
slightly  different  formulation  in  reference  24.  He  points  out  there  that 
n 

if  A  =  H  .a, 

1  '  '  1 
where  is  a  set  of  numbers. 


then  (from  the  Schwarz  inequality) 
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Relation  between  Eulerlan  and  Lagrangian  Microscales 


Taylor  (reference  2)  inferred  an  approximate  relation  "between  Xf] 
and  X  'by  neglecting  the  effect  of  viscosity  on  pres sure -gradient 

|l/2 


3p 


^  W/ j 


w) 


as  "bein^  approximately 


fluctuations  and  estimating 

1/2 

This  led  to  a  constant  ratio  for  all  turbulence. 


The  rough  nature  of  this  analysis  induced  Heisenberg  (reference  25)  to 
conduct  a  more  detailed  study  of  the  static-pressure  fluctuations  and  to 
reestimate  the  X^/X  ratio.  However^  he  followed  Taylor  in  ignoring 

the  dv^d-q-term  in  the  relation  between  and  (see  equa¬ 

tion  (23))  and  in  neglecting  viscous  terms  in  the  relation  between 


Although  these  omissions  are  probably  not  serious  except  in  the 
low  Reynolds  number  range,  it  seems  interesting,  if  only  for  the  sake 

of  completeness,  to  use  a' Heisenberg  type  of  approximation  for 

and  to  repeat  his  treatment  with  the  omissions  rectified. 

From  the  complete  Navier -Stokes  equations  in  the  wave-number  space, 

Heisenberg  deduced  an  approximate  expression  for  (vp)  in  terms  of 
mean  quadruple  products  of  the  "harmonics"  of  the  velocity  field.  His 
principal  simplifying  assumptions  were: 

(a)  Different  Fourier  components  of  the  velocity  field  are 
uncorrelated 


(b)  The  turbulent  energy  spectrum  is  given  by  the  solution  to  his 
equilibrium-energy-transfer  equation,  above  a  lower  cut-off  wave  num¬ 
ber  k^ 

Following  these,  but  using  Chandrasekharas  (reference  26)  solution 
to  the  Heisenberg  equation  instead  of  the  interpolation  formula  used  by 
Heisenberg,  there  results 


where  k  is  a  dimensionless  empirical  constant. 


(27) 
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The  numericai  constant  in  equation  (27)  would  perhaps  have  been 
given  more  accurately  by  the  use  of  a  "self -preserving"  spectrum 
(calculated  by  Chandrasekhar  from  Heisenberg's  equation)  instead  of  the 
stationary  spectrtim  with  low  cut-off  wave  number.  Time  was  not  taken 
to  make  the  requisite  additional  calculations  because:  (a)  'The  value 

of  depends  principally  upon  the  high-wave -number  region  of  the 

velocity  spectrum  rather  than  the  low -wave -number  region,  where  the 
difference  would  be  greatest,  and  (b)  the  experimental  results  and 
(especially)  the  value  of  k  both  have  a  considerable  range  of 
unc  ertainty . 

The  mean  square  of  the  y-component  of  the  Navler -Stokes  equation 
will  lead  to  a  relation  between  and  X: 

dv  Sv  5v  5v 

—  =  —  +  u  —  +  v  —  +  w  — 

dt  ot  ox  oy  oz 


1 

P 


dx 


ay 


dz 


I 


2 

+  W  V 


therefore^ 


(if '  ^(1)'  - 

where  the  correlation  between  pressure  gradient  and  velocity  Laplacian 
function  is  zero  because  of  isotropy. 

Equations  (23)  and  (27)  give  the  first  two  terras  in  terras  of  the 
microscales^  and  the  mean- square  Laplacian  function  is  expressible  in 
terms  of  the  fourth  derivative  of  the  Von  Karman-Howarth  f(r)  correla¬ 
tion  coefficient  at  r  =  0  (reference  27): 


(V2v)2  =  ^(^,)2fiv(o)  (29) 

Consequently,  equation  (28)  becomes 

^  -  (iff  =  ^  ^  ^ "  f 
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The  second  term  in  this  equation  can  he  replaced  hy  the  turbulence 
decay  equation: 


whence: 


iizpl  =  .lov 


=  -5v 


so  that: 


2  +  25v2  =  2|i2  ^  (v^  ^  ^  v2(v')2fiv(o) 


therefore^ 


12.6  ^  12.5 


.  f  fino)  ^ 


Where  =  — 


Batchelor  and  Townsend  (reference  I9)  have  deduced  an  expression 
for  f^^(O)  which  is  valid  in  the  region  of  decay  where  both  l/(v')^ 
and  X^  increase  linearly  with  t  (corresponding  to  large  R\) : 

xMv(o)  =  ^  +  |rxS  (32) 

^  ^  /'av\3/f7sT^l3/2 

where  S  -  J  ^  skewness  factor.  Their  experimental 

results  showed  S  =  0.39  approximately  constant  for  isotropic  turbulence 
Then  the  estimate  for  X/X-n  becomes 


Bx^  I  " 


K 


(33) 
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The  value  of  k  was  first  estimated  by  Heisenberg  (reference  25)^ 
from  measurements  of  turbulence  decay,  as  O.85.  This  method  may  be 
regarded  as  emphasizing  the  (relatively  low  wave  number)  energy -bearing 
range  of  the  spectrum.  Lee  (reference  28)  worked  out  an  estimate  based 
upon  skewness  factor  (k  0.13)^  which  gives  heavy  weight  to  the  high- 
wave-number  range.  Proudman  (reference  29)  has  reestimated  k  by 
comparison  with  measured  curves  of  the  double  and  triple  velocity  correla¬ 
tions.  The  value  k  =  0.45  leads  to  reasonably  good  agreement  for  the 
moderately  high-wave- number  region,  over  a  wide  range  of  R^. 

It  may  be  remarked  that  the  supposed  constancy  of  k  is  merely  a 
postulate  of  the  Heisenberg  dimensional  formulation  of  the  spectral 
transfer  function.  In  fact  it  is  by  no  means  obvious  that  this  turbulent 
part  of  the  transfer  is  quantitatively  independent  of  the  amount  of 
spectrally  local  dissipation  to  heat.  In  any  case,  Proudman* s  estimate 
of  K  =  0.45  has  been  used  here.  Therefore, 


(34) 


In  the  limit  of  R^ — >0,  equation  (34)  does  not  apply  since 
equation  (32)  does  not  apply.  However,  the  appropriate  limiting  rela¬ 
tion  can  be  obtained  directly.  In  this  limiting  condition  the  pressure 
terra  in  equation  (28)  is  negligible  compared  with  the  viscous  terra  (the 
former  goes  with  I/Rx,  “the  latter  with  I/Rx^),  and  the  Eulerian 


velocity  correlation  coefficient  is  (reference  27)  f(r) 
This  gives  f^^(O)  =  3/^^  and 


^-r2/2\2 


(35) 


Eulerian  Analysis  of  Heat  Diffusion  from 
The  two-dimensional  turbulent-heat-transfer 


a  Line  Soiirce 
equation  is 


U 


de 


+  V 


80 

5y 


k 


(i3v) 


(36) 


where  © 
mean,  k 
pressure. 


is  mean  temperature,  -6  is  temperature  fluctuation  about  the 
is  thermal  conductivity,  and  Cp  is  specific  heat  at  constant 
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For  the  thermal  wake  behind  a  line  source  in  isotropic  tiorbulence 
with  constant  mean  velocity  V  =  0.  With  restriction  to  low  turbulence 
levels  a  "boundary -layer"  type  of  approximation  can  be  applied  to  the 
mean  wake,  so  that 

ax2  ay2 

and 


a 

Sx 


(•au)  «  u 

OX 


so  that  equation  (36)  takes  the  approximate  form 


y  5©  ^  k  _  S 

^  ”  P^P  Sy2  "  ^ 


('6v) 


(37) 


It  must  be  emphasized  that  for  this  particular  initial  condition 
on  the  temperature  (effectively  a  "point  source"),  the  restriction  to 
small  turbulence  level  v^/u  «  1  does  not  imply  that  iS*/©  is  small. 
In  fact,  for  this  problem  t9*/0  is  often  greater  than  unity,  especially 
at  the  "edge"  of  the  mean  wake,  as  has  been  discussed  in  reference  10 
and  will  be  brought  out  again  later  in  the  present  report. 

When  the  molecular  transport  can  be  neglected  relative  to  turbulent 
transport 


<38) 

an  equation  given  in  reference  10;  a  slightly  more  general  treatment 
follows : 

With  a  constant  rate  of  heat  generation  (similar  to  steady  state 
in  the  average),  the  application  of  a  Von  Karman  integral -relation 
treatment  to  equation  (38)  yields  an  Integral  condition: 


0  dy  =  Constant  =  H 


(39) 
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where  H  is  the  average  time  rate  at  which  heat  crosses  all  planes 
perpendicular  to  U  per  unit  length  of  heat  _source  (z-direction) .  Of 
course,  ^  has  been  neglected  relative  to  ©U  in  equation  (39) • 

Equation  (38)  has  two  ^mknowns,  and  the  first  objective  is  to 
express  ■5T  as  a  function  of  the  (more  easily  measurable)  ©(x,y). 
After  integration  with  respect  to  y. 


•Qv 


dy  +  F(x) 


But,  by  symmetry,  ■flv  =  0  for  y  =  0.  Hence  F(x)  =0  and 


(i^O) 


This  relation  is  sufficient  for  the  computation  of  ■av(x,y)  from 
the  measured  0(x,y)  but  the  empirical  fact  of  simple  geometrical 
similarity  in  ©(x,y)  suggests  exploitation  of  the  consequent 
simplification. 


Assume 


0(x,y)  =  0o(x)f(|)  \  (4l) 


where  |  =  y/Y'(x).  This  transforms  equation  (39)  to 


0Q(x)y'(x)  = 


H* 

U 


(42) 


where  H*  = 


H 


^00 


=  Constant.  It  transforms  equation  (4o)  to 


2pc^  f(|)  dl 


i 


■I 


^v  =  (-U  I  f(l)  dljY'(x)  ^  +  [u  I  ||  dl)0o(x)  ^  (43) 


dx 
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With,  equation  (^2)^  eliminated  from  equation  (^3)^ 
and  after  integration  of  the  dimensionless  integrals  this  leads  to  the 
final  form  for  the  turbulent -heat -transfer  correlation^ 


^v(x,y)  =  ^  (44) 

where 

I  = 

The  same  sort  of  analysis  can  he  made  on  equation  (37)  which 
includes  the  molecular  conduction^  but  the  rather  large  experimental 
scatter  in  the  present  measurements  seems  to  make  such  a  refinement 
inappropriate - 

An  "exchange"  coefficient  or  "diffusion"  coefficient  for  turbulent 
heat  transfer  kp  is  simply  expressible  in  terms  of  0(x^y).  A  con¬ 
ventional  procedure  for  semiempirical  analyses  is  to  write  for  the 
turbulent  transport  an  expression  just  like  that  for  the  molecular 
transport : 

-pCpiv  =  k;^,  ^  (45) 

which  serves  as  the  definition  of  kp. 

For  the  simple  case  of  equation  (44),  it  turns  out  that 

This  has  the  particularly  interesting  property  that  in  a  nondecaying 
turbulence  at  very  large  values  of  x,  where  the  mean  thermal  wake 
spreads  parabolically  (reference  1),  kp  becomes  independent  of 

explicit  dependence  on  x. 


A  more  startling  simplification  follows  for  a  partic-ular  mean 
temperat*ure  distribution  across  the  wake:  All  dependence  of  kp  on  y 
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disappears  If 


df 


If 


n 


that  is^  if 

is^ 

f(|)  .  e  2  (*•7) 

But  this  is  the  Gaussian  function,  which  is  found  empirically  to  fit 
all  the  measurements  within  the  experimental  scatter.  Hence  one  arrives 
at  the  empirical  result  that  in  both  nondecaying  and  decaying  turbulence 
k/p  is  independent  of  y  in  the  thermal  wake  behind  a  line  source  of 

heat.  From  equation  (47),  M  =  -?f(?)  and 


k^  =  pCpUY'(x) 


dY^ 

dx 


(48) 


It  can  be  seen  that  in  the  nondecaying  case  at  very  large  values  of  x, 
kjji  is  constant  and  independent  of  both  y  and  x. 


Relation  between  Some  Lagrangian  and  Eulerian 
Parameters  in  Transport 

There  has  apparently  been  little  effort  to  relate  the  Eulerian  and 
Lagrangian  formulations  of  turbulent  diffusion  up  to  the  present  time. 
Exceedingly  simple  boundary  conditions  permit  some  connection  to  be 
made  in  restricted  ranges  of  the  present  problem. 

For  nondecaying  or  (with  much  less  accuracy)  decaying  turbulence, 
equation  (13)  applies: 


y2 


=  2 


(n  - 


or,  equation  (12), 
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2  ' 


For  small  values  of  ri.  R_  =  1  -  -li_  whence 

^  X.  2 


y2  =  n^/i  .  1  Jli 


^  =  2n 

dT) 


) 

('■  w) 


(i^9) 


(50) 


SulDstituted  into  equation  (44)  these  give 


^v  =  H*f(|)^^l 


1  jr\v[ 

12  X^2ju 


(51) 


while  equation  (46)  becomes 


krp  — 


pcpv'y  t^(l  -  I  ^ 
dl 


(52) 


where  t) 


In  fact^  for 


small  values  of  t],  ti  —  x. 

U 


For  low-level  nondecaying  turbulence,  ti  =  —  x,  and  equation  (51) 
becomes  TJ 


dv  =  H*f( 


1  v2  x2 


U 


12  u2 


while  equation  (52)  becomes 


(53) 


kip  — 


pCpV 


.y  miA  .  1 1 


v^  x^ 


^  Xr,2 

A.T1  ^ 


(54) 
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As  T) — >■  0  both  equations  (51)  and  (53)  reduce  to 

i  X 

V 

while  equations  (52)  and  (5*^-)  reduce  to 


■Sv  =  H*f(|)-p- 


km  =  DC  V '  ^  ^  V 

W 

At  the  other  extreme,  when  t)  is  very  large, 


where  Mj  = 


J 


=  2tiL^  -  2Mi 
t]Rt^(t])  drj  =  Constant  and 


In  this  case^ 


and 


dY^ 

dT\  ^ 


dv  =  H*f(0 


/ 

0  3/2 

2^T1Lti  -  M; 

kr  -  pkpV 

d| 


ylTj 


2(tiI^  -  Ml) 


1/2 


As  indicated  previously^  the  7]-variation  can  he  expressed  in 
of  X  or  t . 


(55) 


(56) 


(57) 


(58) 


(59) 


(60) 

terms 
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If  T]  is  allowed  to  become  large  enough  to  make  »  M]_ 


dv  =  H*f( I) 


V' 


23/2,3/2^1/2 

kj  =  pcpv-  ^ 


(61) 


(62) 


d| 

The  above  formulas  take  on  particularly  simple  forms  if  the 
empirical  result  of  a  Gaussian  f ( |)  is  utilized: 


f(l)  =  e 


1^2 


Then  the  general  expressions  for  and  krp  (equations  (44)  and  (46)) 

become 


dv(x,y) 


H*  dY 


(Y')2(x)  dx 


ye 


-I/jlV 

2VyV 


V' 

U 


(63) 


and 


kij.(x^y)  =  pc  v'Y'(x) 


dY' 

dx 


(64) 


the  latter  having  been  deduced  in  the  previous  section.  The  particular 
forms  for  small  values  of  r\  would  follow  from  substitution  of  equa¬ 
tions  (49)  and  (50)  into  these  two. 

However^  the  most  interesting  form  occurs  for  very  large  values 
of  r\.  There  Y*(x)  is  given  by  equation  (57) 


i3v(x,y) 


H®^Ly  ^ 
U 


n  37i^ 


If  y 


(tjL,  -  Mj) 


(65) 


=  pCpV^Lp^  -  Constant 


(66) 
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For  still  larger  values  of 


^v(x,y)  = 


and 

kii.(x^y) 

The  constancy  of  kip  for  large  values  of  r|  (large  values  of  t 
or  x)  is  to  he  expected}  a  treatment  of  molecular  diffusion  hy  this 
method  must  certainly  yield  a  constant  coefficient  for  times  much  larger 
than  the  mean  free  time  (of  flight)  of  the  molecules  -  that  is^  for  all 
"macroscopic  times."  Put  another  way^  the  simple  parabolic jDehavior 
of  Y*  for  large  values  of  T]  is  a  sure  indication  that  ©  obeys  the 
simple  classical  diffusion  equation  with  constant  coefficient^  when 
viewed  extremely  "coarsely." 

Perhaps  the  chief  interest  of  equation  (66)  is  its  identification 
of  as  a  significant  Lagrangian  length  for  diffusion  at  a  large 

distance  from  the  source.  It  enters  the  expression  for  turbulent- 
diffusion  coefficient  in  much  the  same  way  as  mean  free  path  enters 
the  expressions  for  the  molecular -diffusion  coefficients.  Furthermore, 
its  role  appears  to  be  much  like  that  attributed  to  Prandtl’s  "mixing 
length,"  which  was  brought  into  the  tiirbulent -transport  problem  in  a 
more  or  less  intuitive  fashion. 

Of  course,  the  possible  crude  nature  of  Taylor’s  original 
T] -postulate  may  render  the  significance  of  Lrj  more  qualitative  than 
quantitative  in  the  case  of  decaying  turbulence. 


T],  such  that  M]_  «  v^lr^y 


H* 


n3/%3/2  (-  i 


(67) 


=  pCpv’Lyj  =  Constant 


COMPUTATION  OF  RESULTS 


Although  Taylor’s  assumption  of  the  unique  dependence  of  -tRt-r 
upon  T}  is  not  likely  to  be  accurate  for  collapsing  together  cases 
with  widely  differing  turbulence  decay  rates,  it  does  provide  a  rela¬ 
tively  simple  relation  between  Y*(x)  and  Rt](t]).  Therefore  all  of 
the  mean- thermal-wake  data  were  reduced  on  the  T]-basis. 

In  principle  the  complete  Rr|(T)  cvorve  can  be  obtained  from  Y^(x) 
by  double  differentiation  (equation  (12)): 
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However,  simple  double  differentiation  of  the  squares  of  a  curve 
as  uncertain  as  Y'(x)  seems  almost  hopelessly  indeterminate  -  although 
Taylor  (reference  2)  and  Collis  (reference  7)  have  apparently  followed 
this  procedure.  A  somewhat  more  circumspect  technique  has  been  tried 
here:  The  values  of  and  ^  were  determined  first,  through  certain 

limit  relations  (to  be  described)  .  Then  the  curve  was  determined 

by  double  differentiation,  subject  to  the  restrictions  of  agreement  with 
the  previously  determined  scales.  Thanks  to  rather  poor  detenninacy  of 
values  of  Xt]  and  this  method  is  not  so  much  of  an  improvement  as 

it  might  first  appear. 


Lagrangian  Microscale  Xtj 

If  equation  (I3)  is  restricted  to  very  small  values  of  n  the 
parabolic  approximation  for  Rt]  can  be  introduced: 


therefore 


(69) 


whence 


(70) 


The  computational  procedlire  was  to  plot  against  and 

to  estimate  the  slope  of  the  faired  curve  at  t]  =  0  where  the  curve 


% 
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must  pass  through  unity.  The  abscissa  intercept  of  the  0-tangent  is 

The  actual  points^  faired  curves^  and'  tangents  for  all  cases  are 
presented  in  figures  13  to  22.  Clearly  the  precision  is  poor. 


Lagrangian  Scale 

Consider  equation  (12)  in  the  limit  as  t] — It  immediately 

gives 


and  the  graphical  procedure  based  on  this  is  also  presented  in  figures  13 
to  22.  Some  of  the  asymptotic  slopes  drawn  are  not  the  best  representa¬ 
tion  of  the  experimental  points.  This  is  due  to  the  auxiliary  (assumed) 
restriction  that  R-q  cannot  increase  with  increasing  values  of  r|  as 
long  as  has  not  previously  dropped  below  zero.  The  graphical 

precision  attainable  is  perhaps  a  little  better  here  than  that  for 

but  the  square  root  necessary  to  get  means  that  X^  is  determined 

about  as  well  as  is  L^. 


Lagrangian  Correlation  Function  RT](‘n) 

With  Xy^  and  L«q  determined^  the  initial  (small  /Sx)  and  final 
(large  Ax)  behavior  of  the  curve  Y^(Ax)  is  prescribed.  These  parts 
of  the  curve  were  drawn  on  a  graph  with  the  experimental  points.  Then 
the  fairing  in  of  a  reasonable  central  portion  to  this  mean  Y*(Ax) 
curve  was  a  relatively  simple  matter.  The  Rt^  ciirve  was  then  obtained 
by  double  differentiation. 

The  curves  drawn  for  Y^(Ax)  in  figures  23  to  32  were  determined 
in  the  fashion  described  above^  as  were  the  curves  for  R^^  in  the  same 
figures . 


Eulerian  Microscale  X 

In  view  of  the  approximate  nature  of  the  determination  of  X^^  no 
new  direct  measurements  were  made  of  the  Eulerian  microscale  X.  Instead^ 
X  was  computed  with  the  energy  equation  for  isotropic  turbulence^  from 
the  measurements  of  turbulence  decay: 
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dt 


(72) 


or,  with  the  space-time  transformation. 


=  -10 V 


v2 


“f- 

dx 


(73) 


Eiilerian  Scale  L 

Earlier  investigations  have  shown  that  the  Eulerian  scale  in  a 
grid-produced  turbulence  is  closely  a  linear  function  of  the  mesh  size 
of  the  grid  producing  the  turbulence  (for  a  given  value  of  x  and  grid 
geometry)  and  is  not  significantly  dependent  upon  the  mean  velocity 
(or  grid  Reynolds  number,  provided  it  is  sufficient  to  cause  turbulence). 
Therefore  the  values  of  L  have  been  deduced  from  earlier  measurements 
at  the  California  Institute  of  Technology  (reference  l8)  on  grids  of 
essentially  the  same  geometiy. 

Table  1  summarizes  the  results  for  Lagrangian  and  Eulerian  scales 
and  microscales.  The  results  have  been  grouped  to  show  the  effect  of 
systematic  variation  of  one  parameter  at  a  time.  Some  of  the  results 
are  presented  in  figures  33^  3^^  and  35* 


Instantaneous  Space-Time  Transformation 

The  permissibility  of  an  instantaneous  space-time  transformation 
in  flowing  turbulence. 


8v 


1  Sv 
U  St 


(7^) 
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can  be  estimated  in  accordance  with  equations  (26).  For  equation  (7^) 
to  be  valid,  the  sufficient  requirements  are  those  given  in  equations  (26) 
that: 


Sv  Sv 


u 


1/2 


«  1 


With  the  aid  of  equation  (23)^  the  turbulence  decay  equation,  and 

the  first  of  these  conditions  can 


the  Taylor  relation 
be  written  in  the  form 


(hvY  _  -  2^ 


1/2 


^\  \  «  1 


(75) 


From  the  Schwarz  inequality,  essentially  the  necessity  that  the 
magnitude  of  any  correlation  coefficient  he  less  than  or  equal  to  unity, 


u- 


^  Sv  <  ,  ,  /  V  /  8v  \ ' 


(76) 


where  the  prime  in  this  expression  denotes  root -mean- square  value.  For 
isotropic  turbulence,  equation  (76)  can  be  written 


(V) 


(77) 
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Thus,  the  second  condition  in  equation  (26)  vill  he  satisfied  if 

2-7  ^«1  (78) 

Both  T  and  v'/u  for  the  flows  studied  are  presented  in  table  1. 
It  is  clear  that  for  these  flows  instantaneous  x  and  t  partial 
derivatives  may  be  taken  proportional  with  reasonable  confidence. 


•flv-Correlation 

The  Eulerian  measure  of  transverse  turbulent  heat  transport  is 
compute^from  the  mean  temperature  distribution.  The  dimensionless 
form,  ■Ov/©qU,  is  given  for  two  typical  cross  sections  in  figures  36 
and  37- 

The  measurements  of  v’/u  and  of  d'/®o  permit  calciaation  of 
the  correlation  coefficient  this  is  also  given  in 

figirres  36  and  37* 

For  the  convenient  and  reasonably  accurate  assiamption  of  Gaussian 
mean  temperature  distribution,  the  corresponding  turbulent-heat-transfer 
coefficient  kT  follows  from  equation  (64).  It  was  found  to  be  inde¬ 
pendent  of  y,  and  typical  curves  of  kip/k  are  given  in  figure  38. 

The  data  for  kip  at  three  different  speeds  behind  the  1-lnch  grid  are 
roiighly  collapsed  together  through  division  of  kp  by  pcpv'l^,  as 

suggested  by  equation  (66a),  an  asymptotic  result  for  nondecaying 
turbulence  (fig.  39). 


DISCUSSION 
Lagrangian  Variables 


Even  a  cursory  examination  of  the  technique  used  in  this  investi¬ 
gation  for  the  determination  of  Lagrangian  correlation  shows  that,  as 
physical  measurements  go,  this  method  is  a  "bad”  one,  largely  because 
of  the  inherent  double  differentiation  between  measured  variable  and 
desired  Information. 

Figures  I3  to  22  suggest  an  uncertainty  in  values  of 
as  large  as  120  percent,  in  spite  of  moderately  good  precision  in  the 
measiirement  of  individual  temperature  distributions  such  as  the  lower 
curve  in  figure  3. 


ho 


MCA  TN  2710 


As  mentioned  earlier^  the  values  of  v’/U  computed  from  initial 
wake  spread  are  consistently  higher  than  those  measured  with  the  hot¬ 
wire  anemometer  (reference  I8) .  The  same  relative  result  was  encountered 
during  a  brief  investigation  following  that  reported  in  reference  I8. 

Up  to  the  present  time  there  has  been  no  satisfactory  explanation  of  the 
discrepancy.  A  tentative  hypothesis  which  would  at  least  account  for 
its  direction  may  be  based  upon  a  human  weakness  in  the  visual  averaging 
of  the  reading  of  a  fluctuation  pointer^  there  seems  to  be  a  tendency  to 
choose  an  "average”  more  or  less  halfway  between  the  extremes  of  the 
needle  travel.  Thus  a  pointer  motion  with  very  skew  probability  density 
(greater  than  O)  would  tend  to  be  "averaged”  at  too  high  a  value.  The 
thermocouple  voltage  in  one  of  these  thermal  wake  traverses  has  just  this 
character  (fig.  11).  Hence  a  visual  averaging  might  yield  too  high  a  wake 
width.  If  this  effect  is  nonnegligible^  it  is  advisable  to  employ  some 
electrical  means  of  averaging  for  skew  signals^  for  example^  the  flux 
meter  and  bucking  circuit  described  in  reference  10. 

In  view  of  the  considerable  uncertainty  in  as  well  as  that 

in  X  and  y'  the  poor  degree  of  agreement  between  experiment  and  theory 
shown  in  figure  33  is  understandable.  Since  the  two  undetermined  con¬ 
stants  in  the  theoretical  result  have  been  evaluated  from  sets  of  experi¬ 
ments  completely  independent  of  the  present  ones^  this  agreement  can  be 
viewed  as  an  affirmative  result. 

Since  some  sort  of  Lagrangian  scale  should  be  a  significant  length 
in  turbulent  heat  and  mass  transport^  as  demonstrated  in  the  analytical 
section  of  this  report^  an  effort  has  been  made  to  find  some  systematic 
variation  in  the  values  of  L^.  Figure  3^  might  be  construed  to  indicate 
a  monotonic  decrease  of  L^/L  with  increasing  values  of  R^.  It  is 
interesting  to  note  that  a  decrease  was  also  observed  for  the  ratio  of 
mixing  length  to  tube  radius  by  Nikuradse  (reference  30)  in  fully  devel¬ 
oped  turbulent  tube  flows.  In  order  to  determine  whether  these  two  rates 
of  decrease  with  increasing  Reynolds  numbers  are  of  the  same  order  of 
magnitude^  an  estimate  has  been  made  of  the  magnitudes  of  Rl  corre¬ 
sponding  to  Nikuradse ’s  results  given  in  figures  28  and  29  of  refer¬ 
ence  30.  Both  scale-to-diameter  ratio  and  average  turbulent  levels  for 
various  Reynolds  numbers  were  estimated  with  the  help  of  Laufer*s  data 
on  tiorbulent  channel  flow  (reference  31)  at  various  Reynolds  numbers. 

The  absolute  level  (i.e.^  the  ordinate  scale)  of  the  resulting  ^max/L 
against  Rj^  c;irve  was  adjusted  to  give  the  most  reasonable-looking  fit 
with  the  L-q/L  data.  This  is  the  dashed  line  in  figure  3^  and  it  shows 
at  least  a  qualitative  resemblance.  It  is  likely^  however^  that  L-irj/L 
is  not  a  unique  function  of  R^. 

Of  course  not  all  of  the  scatter  in  figures  3^  and  35  (involving  Lt|) 
can  be  attributed  to  simple  lack  of  experimental  precision.  Some  is 
evidence  of  the  fact  that  the  Taylor  postulate  of  Lagrangian  correlation 
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function  "being  uniquely  a  function  of  t]  is  certainly  not  very  closely 
true.  Furthermore^  table  1  does  show  rather  systematic  variations  of 
in  some  of  the  three-point  groups.  Most  noticeably^  there  is  a 

regular  decrease  in  with  increasing  Xq/M  (or  perhaps  with 

decreasing  v^U)  ior  each  of  the  three  grids. 


Temperature-Fluctuation  Field 

Fairly  close  behind  a  line  heat  source  in  turbulent  flow^  the 
random  pulse  nature  of  the  temperature  fluctuations  at  a  fixed  point 
has  been  established  by  the  oscillogram  in  reference  10.  This  is 
confirmed  by  the  first  two  oscillograms  in  figure  11^  with  their  highly  . 
skew  probability  densities  at  Ax/m  =  10.  One  of  the  objectives  of  the 
present  investigation  was  to  find  out  whether  this  distinctly  pulsed 
character  persisted  far  downstream  or  whether  molecular  heat  conduction 
becomes  increasingly  effective  in  smearing  out  the  pulses^,  until  they 
are  no  longer  distinguishable  as  such.  The  third  oscillogram  and  proba- 
density  in  figure  11  (Ax/M  =  70)  does  show  a  decided  trend  away 
from  the  pulse-type  signal.  The  molecular  broadening  of  the  laminar  wake 
(corresponding  to  the  pulses)  decreases  the  relative  spacing  of  the  pulses 
in  T9(t)  at  any  point  in  the  "turbulent  wake"  region.  This  is  a  reduc¬ 
tion  in  relative  length  of  the  flat  ('3  =  O)  base  lines  between  pulses 
giving  greater  statistical  symmetry  in  '^(t)  about  its  mean^  that  is^ 
reducing  the  skewness  of  P(‘6). 

A  simple  analysis  will  show  the  existence  of  an  asymptotic  behavior 
of  molecular-conduction  effects  in  a . nondecaying  turbulence.  For  a 
nondecaying  flow  turbulence  and  very  large  values  of  t(~  x/u)  the  mean- 
square  wake  spread  due  to  turbulent  motion  is 


y2s;2Lt|^x  (79) 

u 

On  the  other  hand^  an  approach  to  molecular  diffusion  through 
Taylor^s  concept  of  "continuous  movements"  gives^  for  any  macroscopic 
distance  downstream,  the  mean- square  thermal  wake  width, 


(80) 

u  ^ 

where  A  is  the  mean  free  path  and  c  is  the  root -mean- square  molecular 
velocity. 


k2 
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From  equations  (79)  and  (80) 


(81) 


For  a  typical  case^  take  L  =  1  centimeter,  v’  =10  centimeters  per 

second,  A  =  6  X  10"^  centimeter,  and  c  =  5  X  10^  centimeters  per 

Y  ’ 

second.  T^en  ~  0-17-  For  people  accustomed  to  thinking  of  molecular 

transport  as  negligibly  small  in  turbulent  flow  (e.g.,  in  shear  flow), 
this  ratio  will  appear  quite  large.  The  values  of  krp/k  plotted  in 
figure  38  also  show  that  at  these  low  turbulence  levels  the  molecular 
thermal  conductivity  is  not  necessarily  negligible  compared  with  the 
turbulent  transport . 

The  temperature-fluctuation-level  distributions  across  the 

wake  (figs.  10(a)  and  10(b))  show  the  same  character  as  that  measured  at 
much  higher  turbulence  level  in  a  jet  (reference  10),  with  somewhat 
lower  minimm  values,  which  are  attributable  to  the  lower  turbulence 
level,  A  rough  evaluation  of  the  behavior  of  the  statistical  variables 
in  this  turbulent  thermal  wake  is  obtainable  by  recalling  that  it  consists 
of  a  randomly  "waving”  laminar  thermal  wake.  If  T3(t)  is  crudely 
represented  by  a  randomly  spaced  sequence  of  identical  rectangular  pulses 
with  height  t3q,  width  j,  and  average  spacing  s,  it  is  easily  seen  that 


This  permits  t3*/0  "to  vary  between  0  and  00  as  s/j  travels  the 
permissible  range  from  1  to  00.  Since  points  nearer  the  edge  of  the 
turbulent  thermal  wake  have  higher  values  of  s/j,  the  behavior  of 
equation  (83)  is  consistent  with  the  experimental  distribution.  If 
the  analysis  were  repeated  with  triangular  pulses,  for  example,  the 
quantitative  estimate  would  doubtless  be  more  realistic.  The  higher 
values  of  ('<5*/^)mln  encountered  at  higher  values  of  v^/u  are  an 
indication  that  for  a  given  width  of  laminar  thermal  wake,  the  higher 
value  of  v'/u  leads  to  a  higher  minimum  value  of  s/j. 
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The  closely  Gaussian  shape  of  0/©o  against  y  has  already  been 
pointed  out.  If 


©o 


=  exp 


- 1 
"(xf 


be  introduced,  there  results 


f. 

—  exp 
®o 

1 

-  1> 

J 

1/2 


(84) 


and 


•^  =  exp 


_  izi 


{Xf 


^0 

—  exp 
0 


1/2 


(85) 


Both  of  these  expressions  have  behavior  consistent  with  the  experiments. 

_ _ form  of  dimensionless  transverse  turbulent-heat-transfer  rate 

t3v/©qU  can  be  deduced  for  small  values  of  Ax  (such  that  R-q  ~  l)  with 
this  pulse  representation  of  :3(t).  In  this  picture,  ^  is  the  correla¬ 
tion  between  a  continuous  random  variable  v  and  a  random  pulse  signal 
*6  which  "fires"  every  time  the  continuous  variable  passes  through  a 
specific  value 


Ax 


^  U 


Therefore 


s  Ax 


U  _ 


(86) 


1+1+ 
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where  s  and  ©  are  functions  of  y. 


With  Gaussian  ©(y)^ 


■Qv  _  y_ 


(87) 


The  direct  comparison  "between  this  crude  picture  and  the  experimental 
results  will  be  confined  to  the  correlation  coefficient  R^y  =  t9v/'3*v. 
This  is  of  particular  interest  in  view  of  the  surprisingly  high  experi¬ 
mental  values.  With  equation  (85)  and  the  fact  that  v’/u  =  Y' 
there  results 


—  — 

11 

> 

—  exp 
0o 

i  iL 

- 1> 

—  — 

-1/2 


(88) 


This  contains  the  undetermined  constant  '>9o/®o^  which  can  be  obtained 
from  any  one  of  several  experimental  results.  Figure  36  includes  one 
plot  of  equation  (88)  with  *190/00  determined  from  equation  (8^^-)  and  the 
experimental  value  of  (^  */®^)iiiin^  plot  with  "^q/b  determined  by 

matching  equation  (88)  with  the  experimental  result  at  y/Y'  =  1.4. 


It  is  also  surprising  to  find  that  the  experimental  (^‘av)inax 

large  values  of  Ax  is  even  larger  than  that  at  small  values  of  Ax. 

This  may  be  due  to  a  considerable  experimental  errorj  the  resistance- 
thermometer  voltage  signal  is  much  lower  here.  Unfortunately  no  relation 
corresponding  to  equation  (88)  has  been  deduced  for  large  values  of  Ax^ 
where  R-q  is  essentially  zero. 


The  criteria  for  Taylor’s  hypothesis  of  the  interchangeability  of 
instantaneous  space  U  and  time  derivatives  (assumed  by  him  to 

depend  only  upon  turbulence  level)  have  been  expressed  in  equations  (75) 
and  (78)  as  functions  of  turbulence  level  R\  and  If  X/Xj^  is 

replaced  by  its  theoretical  expression  (equation  (34))  in  terms  of  R-q, 
equation  (75)  becomes 


«  1 


(89) 


For  estimates  of  most  flows  the  first  term  in  the  parentheses  can  be 
neglected)  values  of  R\  less  than  5  or  10  are  rare. 


1^5 
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Since  equation  has  now  been  roughly  verified  hy  experiment, 

equations  (89)  and  (78)  may  serve  as  approximate  criteria  for  the 
validity  of  Taylor's  hypothesis. 

In  the  limit  of  — ^0  when  equation  (35)  replaces  equation  (3^)^ 

there  follows  a  simpler  criterion  to  replace  equation  (75): 

«  1  (90) 

RxVuy 


SUMMARY  OF  RESULTS 


The  following  results  were  obtained  from  the  investigation  of  the 
diffusion  of  heat  from  a  line  source  in  isotropic  turbulence. 

1.  The  thermal  wake  behind  a  heated  wire  set  perpendicular  to  a 
flowing  isotropic  turbulence  (at  sufficiently  low  wire  Reynolds  number) 
consists  of  a  randomly  "waving"  thin,  laminar,  thermal  wake  whose  varia¬ 
tions  in  lateral  position  give  what  may  be  called  the  turbulent  thermal 
wake.  At  a  fixed  point  not  too  far  behind  the  wire  d(t)  is  a  random 
pulse  function,  and  the  nature  of  the  turbulent  heat  transfer  can  be 
deduced  on  this  basis.  Farther  downstream  the  distinct  pulse  nature 
tends  to  disappear. 

2.  The  mean  transverse  temperature  distribution  ^(y)  appears  to 

”0 

be  Gaussian  within  the  experimental  precision  for  all  distances  behind 
the  wire. 

3.  An  Eulerian  analysis  of  this  turbulent -heat -transfer  problem 
permits  computation  of  the  turbulent-heat-transfer  coefficient  k^ 
which  is  essentially  constant  with  respect  to  the  distance  in  the 
direction  of  the  measured  diffusion  y  for  these  boundary  conditions. 

It  is  found  that  at  low  tiu-bulence  levels  (approximately  equal  to  1 
to  2  percent)  the  molecular  heat  transport  is  not  vanishingly  small 
compared  with  the  turbulent  heat  transport. 

4.  Although  Taylor's  postulate  that  Lagrangian  correlations  in 
decaying  turbulence  can  be  made  similar  by  introduction  of  an  independent 

pt 

variable  =  v'(t)  dt  (where  t  is  time  and  v'  is  the  root— 

mean-square  instantaneous  velocity  fluctuation  in  the  y-dlrection) 
seems  to  be  an  oversimplification,  it  has  been  applied  here  for 
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convenience  in  the  reduction  of  data.  .  A  simple  comparison  of  Eulerian 
and  Lagrangian  analyses  for  diffusion  in  nondecaying  turbulence  shows 
that  for  large  values  of  the  distance  from  the  heat  source  Ax  the 
Lagrangian  scale  1^1  enters  the  expression  for  lop,  the  turbulent- 
heat-transfer  coefficient,  much  like  the  empirical  mixing  length  in 
the  old  turbulent  transport  theories.  Therefore  some  properly  modified 
generalization  of  Taylor’s  ii-postulate  should  prove  useful. 

5-  A  correction  and  generalization  of  Heisenberg's  theoretical 
expression  for  the  ratio  of  Eulerian  to  Lagrangian  microscale  x/X-q 
as  a  function  only  of  the  turbulence  Reynolds  number  based  on  microscale 
has  been  made  and  seems  to  agree  roughly  with  experiment.  It  must 
be  noted  that  since  depends  only  upon  a  transformation  dt)  =  v'  dt, 

and  not  upon  the  integral  postvilate  stated  above,  its  validity  is  not 
Impaired  by  any  failure  of  the  integral  postulate. 

6.  Taylor's  hypothesis  for  the  Interchangeability  of  space  and 
time  derivatives  at  low  turbulence  levels  has  been  expressed  in  terms 
of  criteria  which  depend  upon  turbulence  level,  Reynolds  number,  and 
X/X'q.  Applied  to  the  flows  studied  here  it  shows  that  in  these  cases 
such  a  transformation  is  permissible.  By  substitution  of  the  theoret¬ 
ical  expression  for  a  slightly  simpler  and  rougher  criterion 

is  derived,  depending  only  upon  txirbulence  level  and  Rx. 


The  Johns  Hopkins  University 

Baltimore,  Md.,  June  5,  I951 
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(b)  Grid-produced  turbulence. 


Figure  2.-  Shadowgraph  time  e:xposure  of  wire  wake. 
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(a)  No  grid  in  tunnel. 
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Figure  4.-  Spread  of  heat  from  a  line  source.  —  =  43.4^  M  =  1  inch. 
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Figiore  9.-  Ten^jerature  fluctmtlons  at  Xq  =  43  inches,  M  =  1  inch,  and 
U  =  25.6  feet  per  second  along  wake  axis. 
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temperature  fluctuations.  Film  speed 


HP- 


2710 


6? 


(a)  =  0.25. 


(b)  =  0.34. 

X 

Figxare  I6.-  Determination  of  and  -^  =  43.4,  M  =  i  inch,  and 
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(b)  Correlation  function 


Figure  24.-  Spread  of  beat  from  a  line  source  and  correlation  fmction  R 
for  M  =  1  inch,  U  =  25.6  feet  per  second,  and  —  =  43.4. 
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X 

for  M  =  1  inch,  U  =  38  feet  per  second,  and  —  =  43.4, 

M 


(b)  Correlation  function  R^. 

Figure  30.-  Spread  of  heat  from  a  line  source  and  correlation  function  R 
for  M  =  i  Inch^  U  =  25. 6  feet  per  second,  and  —  =  86.8. 
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Figure  33.-  Ratio  of  Lagrangian  and  Eulerian  microscales. 
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Figure  35«-  Hatio  of  Lagrangian  microscale  and  scale. 


(^)  against  y/X*.  Dashed  curves  are  resiilt  of  analysis  based  on 

rectangular  temperature  pulses. 

Figure  3^."  Heat-transfer  correlation  across_thernial  wake,  coii5)uted  from 
measured  mean  temperature  distribution.  U  =  25.6  feet  per  second,  - 
M  =  1  inch,  Xo  =  43.4  inches,  and  Ax  =  10.5  inches. 
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